A terahertz molecular switch 
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We present time-dependent results describing the current through a molecular device, modeled 
as a complex with two active centers connected to leads under bias. We show that, at a properly 
adjusted external voltage, a passing terahertz electromagnetic pulse may cause a transition between 
states of finite and negligible current, suggesting that the system might be useful as a nanoscopic 
switch in the terahertz range. A phase diagram defining the bias region in which the transition takes 
place within a short time is given. As described, the physical processes involved are of an entirely 
different nature than those in ordinary photodetectors. 

PACS numbers: 73.40.Gk, 73.63.-b, 85.65.+h, 74.78.Na 



Individual molecules as well as molecular complexes 
are increasingly being perceived as possible electronic de- 
vice elements, constituting great promise in the process 
of miniaturization. Some successes have already been 
reported along this linei^*^. Among recent examples 
are molecules that exhibit big on-ofF current ratios and 
a large negative differential resistance in two-terminal 
transport^'®, behaving much in the way some mesoscopic 
semiconductor tunneling heterostructures do^^. Such 
performance may be understood decomposing the sys- 
tem into coupled molecular subunits that transit collec- 
tively through a conducting quantum resonance as the 
bias along the device is modified2ii2-ii. The current- 
voltage profile is then determined by the properties of 
this resonance, which in turn may depend strongly on the 
dynamic accumulation of electronic charge in the various 
subunits, thus affecting the tunneling rate for transport. 

When a pair of neighboring subunits of a different kind 
are involved, an energy level associated with one site 
must reach alignment with one in the other site in or- 
der for a device resonance to be established. The narrow 
voltage range in which significant current flows in the ex- 
perimental data of Chen et al^ suggests that such mech- 
anism may be at work in some molecules. Using this two 
site model, Han et al found qualitative agreement be- 
tween the model predictions and the experimental data™. 
A time dependent approach for the analysis of the same 
model allowed us to fit the data to an impressive de- 
gree of accuracy, and showed in addition that switching 
between distinct current states is possible, as reported 
below. Switching is extremely fast compared to that in 
organic complexes, where the effect has been assigned to 
local structural disorder induced by the external electric 
fieldi^. While in the latter switching times are of the or- 
der of 0.1 ^s, we here report characteristic times in the 
order of 0.1 ps, six orders of magnitude shorter. The de- 
vice is then sensitive to terahertz pulses, suggesting that 
the system is potentially useful as a nanoscopic detector 
of electromagnetic radiation in the THz range. 

We consider a molecular complex connected to leads 



subject to an external dc bias potential. As shown 
schematically in Fig. 1, the complex is modeled by two 
active sites, connected to left and right particle reservoirs 
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FIG. 1: Scheme of the device, a, the molecule is placed be- 
tween leads at chemical potentials and fjia- It is charac- 
terized by two active sites with energies Eq, el, which under 
bias, b, change their relative positions allowing for the reso- 
nance condition to be established in a certain bias range. The 
shaded areas show the amount of electron charge at each site. 



at chemical potentials /i^ and /i^, respectively. The de- 
vice is described by the Hamiltonian, 



2 



i=^0,l;a i (#04) -a a=OA;cr a 

£T (T a — 0,1 



(1) 



where {cia ) creates (destroys) an electron with spin 
a at site i and = Ciu represents the occupa- 
tion of site i. The molecular sites are labeled and 1 
and have intrinsic energies Epand £°, respectively, and an 
intra-molecular coupling g. The coupling strength of the 
molecule to the left (right) lead is ^^(ii:), while t repre- 
sents the hoping parameter within the leads. The total 
site energy £i includes the fixed intrinsic energy at the 
site, the external radiation-induced voltage and the ap- 
plied dc bias, the latter represented by a classical term 
linear in the spatial coordinate i. Finally, the effect of 
charge build-up at the molecular sites is accounted for 

I 



by a non-linear term with coupling constant U, the last 
one in the above expression, that leaves out occupation 
with the same quantum numbers. The time-dependent 
wave function that evolves according to this Hamiltonian 
can be expanded in a tight-binding basis as 

|^;c.(T))=5]at(r)|0.,), (2) 



where \(j)ia) is a Wannier state of spin a localized at site i, 
and the coefhcients a*g.(T) obey the non- linear equations 



= £,(r)4,,+i(a^'_i., + 4+i,J (jV -1,0,1,2), 

-1(2), o- t \ k k , k 

- '^-i(2)(,Tja_i(2),CT + ^'L(fl)ao(i),CT + *a-2(3),o-' 
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In this expression the average site densities equal (no) = 

Efc.aho.^l '("1) = Y.k,cr\^i,'y\ ' "^^^^^ ^^'^ sum over 
k,a covers all occupied electron states. In writing Eqs. 
(3) we have adopted a Hartree model for the electron- 
electron interaction. As we will show in what follows, the 
terms proportional to the site densities generated by this 
interaction, nonlinear in nature, are of key importance in 
the behavior of the system. 

Solutions are found discretizing equations (3) and us- 
ing a half-implicit numerical method, which is second- 
order accurate and unitaryi^ii^. It is assumed that the 
coefficients outside the structure are given at incident 
wavenumber k by the expressions 

= {I e'^^" + Rj e-'^^") e"^^^/'' {ja < -L) (4) 

al^=Tj e'^'^'^e-'"'''^ ija>L), (5) 

where the device is explicitly defined to lie between 
sites —L and L, a is the lattice constant, and k' = 



[2m{e — V)/h^]^^^ is the wavenumber of the transmit- 
ted particle, with V = iil ^ fJ-R- The incident amplitude 
/ is supposed to be spatially constant. Also, far from the 
molecule the amplitude of the reflected and transmitted 
waves Rj and Tj are supposed to be weakly dependent 
on site j. This permits to restrict the dependence to the 
linear term, which is found to be an adequate approxi- 
mation provided the time step taken to discretize equa- 
tions (3) is less than certain limit value that depends 
upon the parametrization of the system. In the numer- 
ical procedure the Wannier amplitudes obtained for one 
bias are used as the starting point for the next bias step. 
Once these coefficients are known the current is calcu- 
lated fromi^, 

= f ^ E lHa^>^+i,.)dk (6) 

^ Jo 

We apply the above formalism to a specific situa- 
tion. An excellent fit to the experimental I-V curve at 
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FIG. 2: Current-voltage characteristic in the absence of ex- 
ternal radiation. Parameter values are given in the text, and 
were chosen to fit the experimental data of Chen et aiS at 
T=60K shown as circles. The inset shows the charge in each 
molecule site at different bias values. 



T=60K of Chen et alA is obtained in the absence of ra- 
diation within this model, by choosing the parameters 
eg = 0.52, e° = 0.92, g = 0.008, t = 2.0, U = 0.8 and 
intrinsic level broadenings Tl — Tn = 0.16, all in units 
of eV. Here ^l{r) — '^v'l(r)P{0), with p(0) the density 
of states at the Fermi energy. The energies at the active 
sites are represented by Si = e° — \V + (5i5 sin(27ri^r)]/2, 
( i = 0, 1), where 5E and v are the amplitude and fre- 
quency of the incoming THz radiation, respectively. En- 
ergies are measured with respect to the Fermi energy at 
the emitter contact (/il)- Figure 2 shows our results, to- 
gether with experimental data points from Ref. 5 shown 
as circles. Notice that the current peak is quite narrow, 
and there is a steep fall in the upper critical voltage edge 
14=2.16 V. As the bias is increased and the first site en- 
ergy approaches the Fermi energy charge begins to flow, 
first into site 0, then to site 1 as well as shown in the 
inset and, schematically, in Fig. lb. The effect of this 
charge is to raise the local potential at each subunit, thus 
affecting the transport coefficient of the overall quantum 
mechanical structure. The abrupt fall at Vc is caused by 
a sudden loss of charge at site 1 , and subsequent break of 
the energy level alignment that allows for resonant tun- 
neling through the molecule. 

Figure 3 shows the time evolution of the current flow- 
ing through the system (Fig. 3a) and the charge density 
at sites (full line) and 1 (broken line) (Fig. 3b) when 
electromagnetic radiation of frequency v = IQ THz and 
strength 6E = 1.8 mV is turned on at time =0. The offset 
bias 5V = Vc — V has the value 0.2mV. Before the radia- 
tion enters the sample the system is in resonance, current 
flows and both sites have significant charge. At t=0 the 
external field is turned on and after a time of the order 
of a picosecond, site 1 looses its charge and the current 
drops to a low value. This behavior illustrates the fast 
switching response of the device as a THz pulse passes by. 
The actual transition time is controlled by the internal 
coupling constant, so that in our case T ~ h/g ~ 0.08ps. 
There are small oscillations driven by the external field. 
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FIG. 3: Time evolution of (a) the overall current and (b) the 
charge density in sites and 1 at z/ = 10 THz. Parameters 
used are as in Fig. 2. 



which we have found to become larger in amplitude when 
the signal strength is increased. The fact that the sys- 
tem undergoes a few coherent oscillations before falling 
suggests that the random fluctuations caused by a finite 
temperature are not capable of overcoming the inertia of 
the device. Notice that in the low current state site 1 
has lost its charge. This feature reveals that the system 
is bi-stable very close to Vc- When site 1 is charged, the 
local potential is dynamically adjusted to keep the sys- 
tem in resonance. A disturbance may cause the fiow of 
charge out of the site, however, with the result that the 
resonance condition is lost and current stops. 

Suppose now that the system is in a given state of high 
current, biased slightly below Vc- The transition to a low 
current state is then possible provided the amplitude of 
the external field is large enough to bring the system to 
criticality. Figure 4 shows the region in which the tran- 
sition happens within 2 picoseconds time (labeled YES), 
separated from that region in which the radiation is too 
weak to cause the current to fall that fast, if at all (labeled 
NO). The frequency is = 10 THz and the bias is shown 
in terms of the offset from criticality 6V = Vc — V. A 
rough fit to the boundary is given by SE = 0.014 x SV^^^, 
all quantities in units of V. Our numerical methods al- 
lowed testing frequencies up to order 1/T, and within 
that range we found the boundary to shift to higher val- 
ues as the frequency increases. In the figure, the filled 
circle is a point in the boundary for v = 20 THz while 
the triangle is for v ~ b THz. This is to be expected since 
in the low frequency limit one should have 6E ^ 5V. 

In summary, we have shown that switching between 
different current states in transport through molecules 
is possible and may be triggered by a passing electro- 
magnetic pulse with characteristic times in the terahertz 
range. Depending on the applied external bias, the pas- 
sage of current is turned off by the radiation pulse, thus 
allowing the device to act as a tunable detector of ter- 
ahertz radiation. Switching occurs because of the exis- 
tence of bi-stable states near the I-V characteristic high- 
voltage edge. A similar bi-stable region is known to oc- 



4 



0.012 



0.008 



CD 
to 



0.004 



0.000 




0.000 



0.004 



e51//t 



0.008 



0.012 



cur in transport through GaAlAs double barrier reso- 
nant tunnehng devicesi^4iS4iLi&. A single active energy 
level localized in the quantum well formed between the 
barriers is responsible for the bi-stability in such case. 
Near the bottom of the emitter conduction band a cur- 
rent carrying state through this level, lifted in energy by 
the charge trapped in the well, may coexist with a state of 
no current at the same bias with the well uncharged. In a 
two-site molecular device as the one treated here the ef- 
fect derives from the dynamic coupling between the sites 
within the molecule, and thus has a different physical 
origin. 



FIG. 4; Boundary separating the region at which switching 
takes place within 2 ps of the arrival of an external pulse, la- 
beled YES, from the region in which the well remains charged 
and current still flows beyond that time interval, labeled NO. 
SV measures the external bias relative to the critical edge Vc, 
and 5E is the radiation strength, both in units of the lead 
hopping parameter t. Squares are for incoming radiation of 
frequency = 10 THz. Filled circle and triangle are for fre- 
quencies V = 2Q THz and 5 THz, respectively. 
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